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Abstract. We present steady-state calculations of self- 
similar magnetized accretion discs driving cold, adiabatic, 
non-relativistic jets. For the first time, both the magnetic 
torque due to the jets and a turbulent "viscous" torque are 
taken into account. This latter torque allows a dissipation 
of the accretion power as radiation at the disc surfaces, 
while the former predominantly provides jets with power. 

The parameter space of these structures has been ex- 
plored. It is characterized by four free parameters, namely 
the disc aspect ratio and three MHD turbulence parame- 
ters, related to the anomalous magnetic diffusivities and 
viscosity. It turns out that launching cold jets from thin, 
dissipative discs implies anisotropic turbulent dissipation. 
Jets that asymptotically reach a high Alfvenic Mach num- 
ber are only produced by weakly dissipative discs. 

We obtained general analytical relations between disc 
and jet quantities that must be fulfilled by any steady- 
state model of cold jets, launched from a large radial ex- 
tension of thin discs. We also show that such discs cannot 
have a dominant viscous torque. This is because of the 
chosen geometry, imposing the locus of the Alfven sur- 
face. 

Some observational consequences of these cold mag- 
netized accretion-ejection structures are also briefly dis- 
cussed. 

Key words: Accretion, accretion discs - Magnetohydro- 
dynamics (MHD) - Stars: formation - ISM: jets and out- 
flows - Galaxies: nuclei - Galaxies: jets 



1. Cold jets from Keplerian accretion discs 

Jets of plasma are observed around all young stellar ob- 
jects of low mass, some galactic obje cts an d active galac- 
tic nuclei (see the review of Livio 1997 and references 
therein). These jets share common properties, namely a 
high degree of coUimation and (at least in some objects) 
evidences of interrelations with an underlying accretion 

1991 ). A cat- 



and extragalactic jets rely on the interaction between the 
accretion disc and a large scale magnetic field (for an al- 
ternative view see Lery et al. 199E). The field extracts an- 
gular momentum from the disc, thereby allowing accretion 
and drives the jet. This magnetized jet will be naturally 
self-coUimated provided a large enough current is asymp- 



totically maintained (Chan & Henriksen 198C, Heyvaerts 



Norman 1989) 



Following Blandford & Payne ( |1982D (hereafter BP82) 
there have been numerous studies of magnetized jet dy- 
namics, with prescribed boundary conditions at the disc 
surface. As a result, the conviction grew that jets are in- 
deed magnetized, but the question of their launching re- 
mained open. Indeed, if jets are expected to carry away a 
substantial fraction of the disc angular momentum, then a 
precise treatment of disc-jet interrelations must be done. 
In particular, no standard accretion disc (Shakura & Sun- 
yaev 1973) could be used as a proper boundary condition. 



disc (Hartigan et al. 1995 , Serjeant et al. 
egory of very promising models explaining both galactic 
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Hereafter, we call magnetized accretion-ejection structure 
(MAES), an object where accretion and ejection are in- 
terdependent. 

In a series of papers (Ferreira & Pelletier 1993| , Ferreira 
& Pelletier |199|, Ferreira [l997t hereafter FP93, FP95 and 
F97, respectively), such structures were investigated us- 
ing two simplifying assumptions: (1) jets are cold, i.e. en- 
thalpy plays no role in their energetics; (2) the magnetic 
torque due to the jets is dominant. Only geometrically 
thin accretion discs (Keplerian) were studied, the turbu- 
lent "viscous" torque being simply disregarded. In fact, 
these assumptions are common to all theoretical works 
dealing with the connection between accretion discs and 
jets (e.g. Konigl |1989| , Wardle & Konigl |1994 Li |1995| , Li 
1996|) . 

Thanks to a self-similar formulation allowing to take 
into account all dynamical terms, a smooth transition 
between the resistive accretion disc and the ideal MHD 
jet was achieved. The necessary conditions for steadily 
launching cold jets were thereby obtained. It was shown 
that an equipartition field is required, i.e. a magnetic en- 
ergy density close to the thermal energy density (on the 
disc midplane). Moreover, the vertical gradient of plasma 
pressure is the only force that can gently expells off matter 
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in the resistive region, against both magnetic and grav- 
itational compression. Above the disc (at typically two 
scale- heights), matter flows along magnetic surfaces and 
enthalpy is indeed negligible with respect to other ener- 
gies. By this process, from 1% to 10% only of accreted 
mass can be ejected. Since the transition from accretion 
to ejection occurs at the disc surface, the parameter space 
is highly sensitive to any approximation made on the disc 
vertical equilibrium (see discussion in F97). 

The energetic budget of a cold MAES is quite signifi- 
cantly different from that of a standard disc. Indeed, the 
available mechanical power is shared by radiation released 
at the disc surfaces and an outward MHD Poynting flux 
that powers the jets. But in a thin disc, all the available en- 
ergy is stored as rotation. Therefore, the dominant torque 
dictates the form and amount of dissipation that occurs to 
that energy. While the viscous torque produces only dissi- 
pation, hence a local heating source, the magnetic torque 
produces both another reservoir of energy (outward MHD 
Poynting flux) and a local dissipation (Joule heating). As 
a consequence, a stationary MAES with a dominant mag- 
netic torque is such that almost all the liberated power 
feeds the jets (FP95). 

This quite obvious result has strong observational con- 
sequences. Whenever jets are observed and believed to be 
disc-driven, then one should observe as well a lack of emis- 
sion from the inner parts of the disc. This lack could be 
interpretated as a "hole" (Rutten et al. 1992 ), correspond- 
ing to the radial extension of the magnetized disc (cold 
MAES). However, remember that this is a direct conse- 
quence of our assumption (2). Would it then be possible 
to get both a significant emission from the disc and jets? 

The aim of this paper is to relax the assumption of a 
dominant magnetic torque, keeping however the approxi- 
mation of cold jets (thin or even slim discs, but without 
a hot corona). In Section 2, we write down the set of lo- 
cal MHD equations governing the whole structure, present 
all relevant parameters and provide general analytical re- 
sults on disc-jet interrelations. In the following Section, 
we describe our numerical method to solve the problem 
and display global solutions, from the disc midplane to 
the super- Alfvenic jet regime. We explore the parameter 
space for magnetically-driven jets from resistive and vis- 
cous discs and compare our results to other models. We 
then conclude by summarizing our findings in Section 4. 



2. Cold, non-relativistic MAES 

8.1. General MHD equations 

Our MAES is composed of an accretion disc settled around 
a central mass Ms, (compact object or young star) and 
threaded by a large-scale magnetic field. The presence of 
such a field could be explained by two different phenom- 
ena: advection of interstellar magnetic field and/or local 
magnetic field generation by a disc dynamo (Pudritz 1981 , 



Khanna & Camenzind 1994, Torkelsson & Brandenburg 



1994, Riidiger et al. 1995). Although a quadrupolar topol- 



ogy often arises from these studies, it must be kept in 
mind that (i) they are kinematic, hence not taking into 
account the magnetic feedback on the plasma motion (see 
Yoshizawa & Yokoi 1993| ) and (ii) neglect a possible pri- 
mordial field that could have been advected along with the 
fiow. Thus, it is quite difficult to infer from these works 
what would be the final magnetic topology in a realistic 
accretion disc. Moreover, it has been shown that under re- 
stricted conditions, a quadrupolar topology could produce 
jets from a Keplerian accretion disc, although much less 
powerful than those from a bipolar topology (see F97, Ap- 
pendix A). We will thus assume here a bipolar topology 
for the disc magnetic field. 

The whole MAES is described in a non-relativistic 
framework, so there is a limitation of our calculations 
when applied to the central parts of an accretion disc 
around compact objets. However, we are primarily inter- 
ested in the interrelations between the disc and its jets. 
Thus, we only need to verify that matter remains sub (or 
even midly) relativistic close to the disc (i.e., until the 
Alfven point). We neglect the disc self-gravity with re- 
spect to the gravitational field produced by the central 
mass and look for axisymmetric, stationnary solutions. 

Thanks to axisymmetry, the vectorial quantities ex- 
pressed in cylindrical coordinates (r, 0, z), can be de- 
composed into poloidal and toroidal components, e.g. 
u = Up + rire^ and B = Bp + Bc^e^. A bipolar config- 
uration allows us to describe the poloidal magnetic field 
as 



B„ — -Va X Grf 



(1) 



where a(r, z) is an even function of z and a = constant 
describes a surface of constant magnetic flux (a — rA^, 
being the toroidal component of the potential vector). 
The following set of equations describes a non-relativistic 
MAES (FP93): 



Mass conservation 
V.pu = 

Momentum conservation 

pu.Vu = -VP - pV$G + J X B + V.T 

Ohm's law and toroidal field induction 

VrnJ^e^ = Up X Bp 
1 

Up — HpQr) 



V.(%VrB^) = v.- 



(2) 

(3) 

(4) 
(5) 



where p is the density of matter, P the thermal pressure, 
<i>G = —GM^,/{r'^ + z^)^/^ the gravitational potential, 
J = V X B//io the curr ent, T the "viscous" stress ten- 
sor (Shakura & Sunyaev 1973 ), = ?7rn/Mo and v!^ the 
(anomalous) poloidal and toroidal magnetic diffusivities. 
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All transport coefRcients appearing in the above set of 
equations, namely Vm, v'^ and ("viscosity", contained 
in T), are assumed to be of turbulent origin. 

Energy conservation is not self-consistently solved here 
(see the discussion in FP95). Since jets are cold, the precise 
way energy is transferred from the disc midplane to the jet 
plays no dynamical role. We can therefore freely specify 
the temperature profile. While previous works assumed 
isothermal magnetic surfaces, we will here use an adiabatic 
equation 



P = K{a)p' 



(6) 



where 7 = f for a monoatomic gas and if (a) (related 
to the specific entropy) is conserved along each field line. 
Finally, in order to close the system, we use the perfect 
gas law 



P = nkuT 



(7) 



where T is the temperature of the plasma, fcs the Boltz- 
mann constant and n = pj irip {nip being the proton mass). 

2.2. Relevant disc parameters 

Taking into account all dynamical effects appearing in 
Eq. (^) to (|^) is costly, for it results in quite a large number 
of dimensionless parameters. However, thanks to various 
constraints, only four of them are really free (but the pa- 
rameter space of a cold MAES remaining to be spanned 
is large, see Sect. 3.3). Below, we list all the relevant pa- 
rameters of a cold MAES. 

The strength of the magnetic field in the disc is mea- 
sured by 



f^oP 



(8) 



This parameter cannot be much smaller than unity for the 
disc would become prone to shearing instabilities (Balbus 
& Hawley 1991). In the present work however, we always 
find II of order unity (rough equipartition) . 
The efficency of ejection, defined as 



din Ma 

dlnr 



(9) 



is the key parameter linking accretion to ejection (mass 
conservation). For ^ = we obviously have a standard 
disc, whereas a typical value for cold MAES is 10^^. This 
parameter is strongly constrained by both the disc ver- 
tical equilibrium (minimum value of ^) and the steady 
production of super- Alfvenic jets (maximum value, F97). 
This local ejection efficency is actually related to the mag- 
netic flux distribution. When the magnetic flux scales as a 
power law of index /3 (i.e., /3 = din a/c? In Tq is a constant), 
the following scaling must be fulfilled in a Keplerian disc: 
P = 3/4 + ^/2 (note that BP82 used /3 = 3/4). 



The (effective) magnetic Reynolds number, defined as 



rUr 



where 
a(r, 2 = 0) 



(10) 



9^ 



z=0 



is a direct measure of the bending of the magnetic field 
lines at the disc midplane (where, from Eq. (^), one gets 
the right handside of the above equation). In previous 
studies, TZm was always of order e~^, hence providing 
a large amount of toroidal current at the disc midplane 
(FP95, F97). This current induces a large radial compo- 
nent of the magnetic field, such that Br Bz at the 
disc surface. This is the required bending for magnetically 
launching cold material (BP82). However, another situa- 
tion with a much smaller TZm (> 2 — /?) could still lead to 
cold ejection, provided there is an extra source of toroidal 
current at the disc surface. 

The ratio of the magnetic torque due to the jet to the 
turbulent "viscous" torque at the disc midplane, namely 



A 



(J X B) 



(V.T), 



(11) 



could, in principle, range from zero (standard disc) to in- 
finity (all previous studies on magnetized accretion disc 
have neglected the viscous torque) . Note that by "viscous" 
torque, we mean a radial transport of angular momentum. 
This "viscous" transport could originate from MHD in- 
stabilities, i. e. fluc tuations of small scale magnetic fields 
(Stone et al. [19961 ). 

The toroidal magnetic field at the disc surface is mea- 
sured by 



h dB, 



Bo dz 



(12) 



z=0 



The shearing of the magnetic configuration is related to 
the current density flowing at the disc midplane and thus 
depends on the global electric circuit. This parameter was 
found to be of order unity for strongly bent poloidal field 
(i.e., Bcf, ^ Br ^ Bz at the disc surface, F97). 

All above parameters are determined by physical con- 
straints. Below, we present the parameters that remain 
free, one geometrical and three related to the MHD tur- 
bulence {vTn,v'r,-^,Vy)- 

(i) The first parameter is the disc aspect ratio 



e = 



h{r) 



(13) 



where h{r) is the local disc vertical (pressure) scale height. 
Having no energy equation, this parameter remains free 
and £ is a constant all over the magnetized disc (FP93). 
With our treatment, we are able to keep all terms in the 
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equations. Therefore, we do not use the usual approxima- 
tion £ ^ 1. In fact, solutions will be obtained for thin 
(e = 0.01) to slim (e = 0.1) discs. 

(ii) The second parameter measures the strength of 
the MHD turbulence 



(14) 



VAh 



where Va is the Alfven speed. If one requires stability 
against resistive instabilities (naively, = h'^ ^ ta — 
h/VA), one gets Um of order unity. 

(in) The third parameter is the magnetic Prandtl mun- 

ber 



which measures the ratio of the "viscosity" to the poloidal 
magnetic diffusivity at the disc midplanc. This parameter 
plays an important role in the angular momentum equa- 
tion, which writes at the disc midplanc 



7'„(1 + A) =7^„ 



(16) 



Our usual understanding about turbulence would impose 
Vm of order unity (Pouquet et al. 1976 ). If this remains 
true in a MAES, there is a strong link between the dom- 
inant torque (A) and the magnetic field curvature (TZm)'- 
strongly bent fields {TZm ~ e~^) extract all disc angular 
momentum. The global energy budget of a cold MAES 
(with e ^ 1) reads 

Pub ~ 2Prad + 2PmHD , (17) 

where Pub is the liberated mechanical power, Prad the disc 
luminosity and Pmhd the outward MHD Poynting flux. 
The liberated power is Pub ^ Pace = GM^.Mae/'^ri, where 
r.i is the disc inner radius and Mae is the accretion rate 
at the disc outer edge Ve (FP95). Using the disc angular 
momentum conservation gives 



IP 



1 



lih 



2P, 



jet 



1 + A 
A 



p.. I + A- 

Thus, a stationary MAES with Vm ^ 1 and a dominant 
magnetic torque is such that only a fraction h/r of the 
total power is released as radiation. However, for the pur- 
pose of a general investigation, we will keep Vm free in 
this paper. 

(iv) The last parameter measures the anisotropy of the 
magnetic diffusivities 



(19) 



z=0 



at the disc midplanc. An isotropic turbulent dissipation 
would be described with Xm of order unity. However, we 
expect the leading instabilities triggered inside the MAES 
to produce enhanced dissipation of the toroidal field, thus 
providing < 1. 



2.3. Launching jets from Keplerian discs 

In order to eject matter from the disc, a necessary condi- 
tion is that the magnetic torque changes its sign at the disc 
surface. Indeed, it brakes the matter in the disc midplane 
and both angular momentum and energy are stored in the 
magnetic field. Above the disc, angular momentum and 
energy must be given back to matter for launching a jet, 
namely (J x B).e0 > at z = /i. The toroidal component 
of the Lorentz force is mainly governed by the behaviour 
of the radial current density, which is in turn controled 
by the induction equation (5). Integrated vertically, this 
equation gives 



(15) 77;„ J, = liJo + r / dzBp.Vr! - B, 



(20) 



The last term (field advection) is negligible in the resis- 
tive disc but will exactly balance the differential rotation 
effect in the ideal MHD region. The change of sign of the 
magnetic torque implies that the radial current density 
must decrease on a disc scale height. After some approx- 
imations]^ the necessary condition for launching a jet is 
found to be 



A = 



'iXr. 



.eV„ 



(21) 



This very important relation does not allow much free- 
dom in thin discs. Indeed, for conventional values of the 
turbulence parameters {um ^ Vm ^ Xm ^ l)i one gets 
A ~ £~^, namely a dominant magnetic torque. To enforce 
comparable torques (A 1), one then must ask for ei- 
ther an anisotropic turbulence (xm < 1) or a large Vm (or 
both). A fortiori, launching cold jets with A ^ 1 seems 
almost impossible. Cold jets have a tremendous influence 
on the disc structure. One useful quantity to evaluate is 
the ratio a of the MHD Poynting flux to the kinetic energy 
flux, 



A 



<7SM 



C(l + A) 



(22) 



measured here right above the disc (at the slow magne- 
tosonic point). It shows that, unless ^ is of order unity or 
A ~ ^, the magnetic field completely dominates matter at 
the disc surface (i.e., asM ^ !)■ 

To summarize, A is constrained by the very existence 
of a jet, TZm by the disc angular momentum conservation 
and q is an explicit functions of the other parameters: 



q 



A 



A 



(23) 



with S = flo/i^K, ratio of the angular velocity at the disc 
midplane to the Keplerian rate fix (see the exact expres- 
sion of S in Eq. (^)). The two remaining disc param- 
eters, ^ and /i, are constrained by the smooth crossing 

^ see Appendix B of F97. In order to obtain a generalization 
of his expressions, replace his condition F ~ 1 by F ~ TZm£. 
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of two critical points (slow magnetosonic and Alfven, see 
FP95). Thus, the present model has only four free param- 
eters. In the next section, we show the strong links be- 
tween cold jet parameters and those describing the disc. 
In practice however, we will impose a set of disc param- 
eters {s,am,'Pm,Xm) and obtain a posteriori the jet pa- 
rameters. 

2.4- Relevant jet parameters 

The jet is the region where the transport coefhcients are 
equal to zero, namely the ideal MHD medium above the 
disc. In this regime, mass and flux conservations combined 
with Ohm's law (0) provide 



77(a) 



(24) 



where 77(a) — y/fioPA is a constant along a magnetic sur- 
face and Pa is the density at the Alfven point, where the 
poloidal velocity reaches the poloidal Alfven velocity. The 
induction equation (ph becomes 



r2*(a) — CI — rj 



' Popr 



(25) 



where r2*(a) is the rotation rate of a magnetic surface 
(very close to the Keplerian value). In the jet, plasma 
flows along a magnetic surface with a total velocity u = 
(7///iop)B -t- ri^re^ not parallel to the total field. 

The angular momentum conservation in the jet simply 
writes 

rB^ 



(26) 



where is the Alfven radius. Above the disc, the turbu- 
lent torque vanishes and only remains a magnetic acceler- 
ating torque. 

The projection of the momentum conservation equa- 
tion along a magnetic surface provides the Bernoulli equa- 
tion 



u 

Y 



rB, 



= £{a) + nlr\ = E(a) 



(27) 



where E{a) is the constant specific energy carried by the 
jet, H is the specific enthalpy defined as H — 7-P/(7 — 
l)p in the adiabatic case. Enthalpy is negligible for "cold" 
plasma ejection at the disc surface. 

The shape of the magnetic surface is given by the jet 
transverse equilibrium, namely the Grad-Shafranov equa- 
tion that can be written in the following form 



(1 - r7z^)j0 = Ja + Jk , 

where the sources of current are 



(28) 



Jx 



pr 
Bl 



d£ 2^-^* 
da da 



dfl^,r\ 
da 



m^B^ din PA 



2 Va 



2po 



da 



m .V Inp 

Por 



Here, — Up/Vj^ is the Alfven Mach number and g = 
l-fl/n^ (Pelletier & Pudritz [1992|, hereafter PP92). 



We choose to use the usual jet parameters as defined 
by BP82. Equation (|2^) allows to define a mass load pa- 
rameter 



K = rj- 



(29) 



where the subscript o defines quantities at the disc mid- 
plane. This dimensionless parameter describes the mass 
fiux per magnetic flux unit {dMj/da = 2^7]/ po), thus it is 
constrained by the disc vertical equilibrium. 

The total specific angular momentum, defined by 



A 



^*7^^ 
5 'o' o 



(30) 



provides a measure of the magnetic lever arm acting on 
the disc, i.e. a measure of the Alfven radius r^. In the 
case where all the MHD Poynting flux is converted into jet 
kinetic power, the magnetic lever arm parameter uniquely 
determines the jet terminal velocity 



r!oro(2A - 3)1/2 



(31) 



The asymptotic jet behaviour is strongly influenced 
by another normalized quantity. This last jet parameter 
( "fastness" parameter) is a direct measure of how fast the 
magnetic rotator is (Michel |1969| , PP92, F97, Lery et al. 



1998) 



UJA 



Vap,a 



(32) 



PP92 showed that super-Alfvenic jets require oja bigger 
than unity. This subtle parameter links rotation of the 
magnetic surface to poloidal motion. In our calculations, 
it is intimately related to the MHD power still available 
at the Alfven surface. For uja bigger than but close to 
unity (or gA < 1/2), matter reaches the Alfven surface 
at the expense of all the MHD Poynting flux (no more 
current available) but recoUimation takes place just after 
it. For larger uja (typically bigger than 1.4, or gA > 1/2), 
a large amount of current is still available and the jet 
propagates much farther away before recoUimation oc- 
curs (F97). The "fastness" parameter lua determines also 
the magnetic geometry at the Alfven surface, both pitch 
(arcsin {—B^/Bp)) and opening {9 — a,rccos{Bz/ Bp)) an- 
gles. The larger loa, the larger these angles (centrifugal 
effect stronger than the hoop stress), providing a subse- 
quent large widening of the jet. The fastness parameter 
UJA is directly proportional to the toroidal field at the disc 
surface q (see below). This parameter did not explicitely 
appear in BP82, but it was hidden in the initial angle of 
field lines at the disc surface (the larger angle, the larger 

In usual treatments of MHD jets, where the disc is a 
mere boundary condition, k is fixed by the regularity con- 
dition at the slow point whereas A is fixed by the Alfvenic 
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condition for a given uja- In our treatment, all these cold 
jet parameters are completely fixed by the set of disc pa- 
rameters. Using ideal MHD equations, mass conservation 
and angular momentum conservation, we get 



A 



A 



2^(1 + A) 



9 , ^XlTlNj^ —1/2 



A- 1 



(33) 



UJA 



q 



A3/2 sm{^A - 6a) 



A- 1 



sm (f>A 



where 4) a is the angle between the Alfven surface and the 
vertical axis, and 9a is the opening angle estimated at the 
Alvenic transition. This angle 9a is implicitly determined 
by the Grad-Shafranov equation (see Appendix A). The 
angle (pA is imposed by the chosen geometry of the Alfven 
surface. Note that this last expression is valid for conical 
Alfven surfaces (see Appendix B), a geometry which arises 
naturally when the above parameters vary slowly across 
the jet. Rewriting the last expression as 



UJA 



A3/2 sin(0A - 9a) 



A- 1 



sm (pA 



(34) 



where av is the usual alph a par ameter for turbulent vis- 
cosity (Shakura & Sunyaev 1973), and using the necessary 
condition for super- Alfvenic jets [oja > 1) provides 



A > 



A-1 1 



A3/2 



(35) 



This is a strong constraint on the underlying accretion 
process. A dominant viscous torque (A < 1) would then 
require either (1) at, smaller than unity but huge magnetic 
lever arms, namely A > £~^, or (2) > e~^. Whether 
the last condition is clearly unphysical, the former remains 
open at this stage. 



3. Global solutions from disc to jets 
3.1. Numerical approach 
3.1.1. The self-similar ansatz 

In order to have a global solution for a MAES, we need 
to solve the whole set of MHD equations (||) to (Q). Be- 
cause these equations are highly non-linear partial differ- 
ential equations, we use a variable separation method that 
greatly simplifies the resolution. Self-similarity allows us 
to write all the quantities like 



Q(r,z)^Q,{-] ' fgix) 



(36) 



where x — z/h{r) and r,, is the outer radius of the mag- 
netized disc. By doing this, we obtain a set of ordinary 
differential equations acting on /q where we can consider 




z(a)/i'„ 



Fig. 1. Characteristic velocities along a magnetic surface, 
normalized to the local disc sound speed ^Kho'- slow 
and fast magnetosonic velocities Vsm and Vfm, critical 
plasma velocity V — u.n (n is the relevant direction of 
wave propagation, FP95) and the Alfven velocity Vau 
(note that V = Vau is equivalent to Up = Vap). These 
curves describe the MAES presented in Fig. ||. However, 
the behaviour is common to all self-similar, cold jets from 
Keplerian discs obtained so far. Note that the jet is super- 
fast in the usual sense {up > Vfm, F97). 

all terms of each equation. The values of the coefBcents 
aq are given by FP93. The magnetic flux is written 



a(r, z) 



ip{x) 



(37) 



To obtain the variation of any quantity Q along a magnetic 
surface, we just have to write 



Qix) = Oo^-"''/^(x)/q(x) 



(38) 



where x must be understood here as a curvilinear ab- 
scissa along the magnetic surface. In fact, we find more 
convenient to choose as a curvilinear variable z{a)/ro = 
exip~^/^ (x), since it directly provides a physical scale once 
To is chosen. 

3.1.2. Turbulent transport coefficients 

In the self-similar framework, only the vertical variation 
of the turbulent transport coefficents must be prescribed. 
Since these anomalous coefficents are expected to arise 
from an MHD turbulence triggered inside the disc, we as- 
sume that the magnetic diffusivities Vm and u'^ decrease 
on a disc scale-height. We therefore used a Gaussian law 
for their vertical profile. 

The same issue appears for the turbulent "viscous" 
torque, which is responsible for the radial transport of 
angular momentum. Following the above arguments, the 
corresponding turbulent viscosity Vy would also decrease 
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on a disc scale-height. However, we must take care to con- 
serve the total disc angular momentum. For this purpose, 
we write the angular momentum conservation equation as 

\7.{pnr^Up-'^Bp-rT) = (39) 

where we used 

1 

-V.(rT) = (V.T).e0 and r = 

MHD discs launching cold jets have comparable radial and 
vertical gradients of the angular velocity S7, hence provid- 
ing comparable components Tr^ and of the turbulent 
stress tensor. Thus, prescribing the vertical profile of the 
poloidal components of r (instead of the torque itself or 
the viscosity), insures the conservation of the total angu- 
lar momentum. Because the turbulent torque must be an 
even function of a;, we use a exp(— a;^) law intead of a sim- 
ple Gaussian law. The resulting turbulent torque remains 
always negative in the disc and decreases on a disc scale 
height. 

3.2. Numerical results 

With a given set of disc parameters {s,am,'Pm,Xm), we 
integrate the set of ordinary differential equations step by 
step with respect to the self-similar variable x — z/h{r). 
The transition between the resistive, viscous disc and the 
ideal MHD medium occurs smoothly above the disc sur- 
face. 

The determinant appearing in the set of ideal MHD 
equations vanishes when the critical plasma velocity is 
equal to the phase velocity of a typical wave of the 
medium: slow magnetosonic (SM), Alfven (A) and fast 
magnetosonic (FM) waves. As already shown by many 
authors (e.g. BP82 and FP95 for the exact self-similar 
expressions of these velocities), the expressions of these 
velocities are modified by self-similarity. The slow point 
is mostly related to the vertical velocity, the Alfven to 
the poloidal one, whereas the fast is related to the radial 
velocity. 

For given (e, a™, Xm), the smooth crossing of these 
points is allowed by fine-tuning two parameters, namely 
jjL and f . For the slow point, the strength of the magnetic 
field jjL must be equal to a critical value jjLc- If /i < /ic the 
density profile is flatter, leading to a fall of the vertical 
speed of matter, whereas if /i > ^c, the magnetic field 
pinches too strongly the disc, forbidding matter to leave 
the disc. The crossing of the Alfven point is controled by 
the efficency of ejection i^. If ^ > the magnetic field 
lines are too much opened because of an overwhelming 
centrifugal effect. This leads to fir^ > ri*r^, that is > 
and the structure is non-steady. On the contrary, £, < 
leads to an unphysical closing of the magnetic surfaces 
due to the magnetic tension. For each attempt to cross 
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Fig. 2. Typical solution from a dissipative disc (A = 2.43) 
with e = 0.1, a,n = 1, V^n = 3.2 and x™ = 0.259. All 
quantities are shown along a magnetic poloidal line, an- 
chored on a radius Tq. The density (upper left pannel) is 
normalized to the density po at the disc midplane and jet 
poloidal velocity (lower right pannel) to the Keplerian ve- 
locity. Stars show the location of the slow-magnetosonic 
(SM) and Alfven (A) critical points. The magnetic sur- 
face (upper right pannel) presents the usual recoUimating 
behaviour once all MHD power has been converted into ki- 
netic power (plateau in the poloidal velocity). The lower 
left pannel shows both radial and vertical components 
of the plasma velocity inside the disc, normalized to the 
sound speed. The transition between the resistive disc and 
ideal MHD jet occurs above the disc surface {z/h ^ 1.6), 
where both components become comparable. 

the Alfven point by fine-tuning ^, we must find again the 
critical fj,. 

As in F97, all trans- Alfvenic solutions presented here 
do not cross the last critical point, the fast magnetosonic 
one. The numerical integration stops at exactly the loca- 
tion of this point (see Fig. |^). The meeting of this last 
point is unavoidable provided the numerical integration is 
done sufficiently far. Actually, it seems inherent to self- 
similar solutions that the crossing of the three critical 
points is impossible (see discussion in F97). 

Only high-LUA jets propagate much farther away than 
the Alfven point, widening a lot {too ^ ta) and reaching 
large asymptotic Alfvenic Mach numbers (m^ » 1). We 
label these jets as being powerful. Such jets are produced 
whenever the initial MHD Poynting flux is high, namely 
when the radial current density is large on the disc mid- 
plane. Thus, a high value of q (close to unity) is required 
for cold jets launched from a large radial extension in the 
disc. Then, such jets exert a large torque on the disc mid- 
plane, thereby producing a large accretion velocity. 

For a small poloidal magnetic diffusivity, the poloidal 
field will tend to be strongly bent by the accretion flow. 
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Fig. 3. Parameter space of cold, non-relativistic magnetized accretion-ejection structures for a,„ = 1 and e ~ 0.1. 
The shaded area corresponds to the location where numerical solutions could be found. Thick dashed lines show 
theoretical limits: super-slow magnetosonic (upper limit) and super- Alfvcnic (lower limit) flow. While the range in Xm 
is quite narrow, the one covered by the magnetic Prandtl number is quite large (from 10^^ to almost 4). Levels for 
the corresponding jet parameters k and loa are also displayed. 



eventually providing a large magnetic Reynolds number. 
In this case, the turbulent torque becomes negligible and 
all the available power goes into the jet. Not much power 
is finally radiated at the disc surfaces. Requiring both 
high-Wyi, cold jets and a dissipative disc is quite a diffi- 
cult task, that can however be achieved in two different 
(though costly) ways. 

Equation (|l^) shows that for a magnetic Prandtl num- 
ber Vm of order unity there is a strong link between the 
magnetic configuration (TZm) and the magnetic torque 
(A). Increasing the magnetic diffusivity only (a™), so 
that the field lines remain almost vertical inside the disc 
{TZm ^ 1), leads to comparable torques. Therefore, a com- 
parable amount of energy is released as radiation and 
MHD Poynting flux. The other possibility consists on in- 
creasing the magnetic Prandtl number only (Vm > !)• In 
this way, the effects of the turbulent viscosity are enhanced 
so that comparable torques can still be achieved, even with 
a strongly bent poloidal magnetic field {TZm ^ s~^). 

Actually, it is noteworthy that using the Shakura- 
Sunyaev prescription for the viscosity {uy = ayflxh^) 
gives 

ay = amVmti^'^ . (40) 



Thus, both ways of obtaining dissipative discs producing 
high-ujA jets imply an increase of the viscosity parameter 

Figure ^ illustrates a typical dissipative solution, de- 
scribed with the set of parameters (e, a^, T'm^ Xm) = (0.1, 
1, 3.2, 0.259). This solution is in the upper right side of 
the parameter space represented in Fig. ^. It corresponds 
to a dissipative disc (A = 2.43, 29% of the total me- 
chanical power being released as radiation) and high-ijjyi 
jets. The resulting set of jet parameters is {k,\,ljja) — 
(5 10-3,78,1.44). 

Inside the disc, both the radial and vertical compo- 
nents of the plasma velocity are negative. Matter is ac- 
creting towards the central object and slightly converging 
toward the disc midplane. The transition from accretion 
to ejection occurs at the disc surface, where all dynamical 
terms are comparable. At the surface, the steep decrease 
of the density profile {p^ = p{h) ~ epo, where po is the 
density at the disc midplane) goes along with a typical 
outflow velocity 

ut ^ Uoi ■ (41) 

This transition occurs in still resistive and viscous upper 
layers, where the magnetic torque azimuthaly accelerates 
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Fig. 4. Parameter space of cold, non-relativistic magnetized accretion-ejection structures for am = 0.3 and 
For a smaller am, the region where solutions can be found is much smaller (see Fig. with a magnetic 
number ranging from 10^^ to 1.2. Here no dissipative solution could be found and jets could not reach an loa 
than 1.15. 
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the plasma. Farther out, both the decrease of the transport 
coefficients and the Lorentz force itself, enforce plasma 
to flow along magnetic surfaces. Once in this ideal MHD 
regime, the flow encounters the first critical point (SM). 

As previously said, the asymptotic jet behaviour is ex- 
actly the same as in F97. Plasma always achieves its max- 
imum velocity (mostly vertical), almost all initial MHD 
Poynting flux being finally converted into kinetic energy 
(Fig. 1^). Jets become super fast-magnetosonic in the con- 
ventionnal sense, namely Up > Vfm — Va<p {B^ ^ Bp). 
From this point on, both the centrifugal force and total 
pressure gradient cannot overcome the hoop stress. This 
leads to a recollimation of the magnetic surface (negative 
opening angle), until the last critical point is finally met 
(Fig. |l|). The maximum radius achieved depends mostly 
on UJA (see Fig. 12 in F97). 

High-Ljyi jets from a dissipative disc were obtained here 
by increasing the magnetic Prandtl number. In this par- 
ticular case, we obtained a phenomenological viscosity pa- 
rameter ay ~ 2.2, bigger than unity. This is problematic 
since one usually expects a turbulence which is both sub- 
sonic and with a correlation length smaller than the disc 
scale height (hence a„ < 1). Alternatively, we did find a 
dissipative solution with a^ = 0.27, but the price was an 
increase of the MHD turbulence parameter am — 1-8. It 



seems therefore too costly to get both high-w^i jets and 
dissipative discs. To firmly settle this issue, we scanned 
the parameter space of cold MAES. 

3.3. Parameter space of adiabatic, cold jets from 
Keplerian discs 

The parameter space of cold MAES is obtained by varying 
the set of free disc parameters (e, am, 'Pm, Xm)- We choose 
to fix the values of both e and am and represent the pa- 
rameter space with the remaining parameters (Fig. ^ and 
. Numerical solutions are only found inside the shaded 
areas, where we also plot levels of two jet parameters 
UJA and K. This region is embedded inside a larger re- 
gion (thick dashed lines), obtained by two analytical con- 
straints. 

The first constraint arises from the requirement that 
jets become super-slow magnetosonic. It is thus related to 
the disc vertical equilibrium: a too strong magnetic pinch- 
ing of the disc forbids plasma to escape from it. Since 
plasma pressure is the only force that allows matter to be 
lifted from the disc, it cannot be much smaller than the 
vertical component of the magnetic pressure. This mag- 
netic pressure is due to the growth of both radial and 
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Fig. 5. Parameter space of cold, adiabatic jets from Keplerian discs with am = 1 and e = 0.1. It illustrates the disc 
parameter space in Fig. The lower limit of the shaded area is imposed by the Alfvenic constraint: no trans- Alfvenic 
jet can be found below. The upper limit is imposed by the disc vertical equilibrium (trans-SM jet). This parameter 
space for jets is to be compared with those obtained by Blandford & Payne (1982), Wardle & Konigl ( 1993| ) and Li 
( |1995| ). Using this plot along with Fig. || allows to compute all other MAES parameters (^, /i, q). 



toroidal magnetic fields. Using a Taylor expansion (on z), 
this constraint can be written as 

Each contribution to the magnetic compression is ex- 
pected to be smaller than the one provoked by the vertical 
magnetic field (equipartition between plasma pressure and 
magnetic pressure at the disc midplane). The upper limit 
appearing in Fig. ^ and ^ corresponds to a necessary con- 
dition involving only the radial component, namely 

Xm < ^(V2-P™e) . (43) 

This Taylor expansion is however not valid in discs where 
the poloidal magnetic field is straight [TZm ^ 1). In that 
case, there is an extra source of toroidal current at the 
disc surface, leading to a large bending of the poloidal 
field. The corresponding magnetic pressure increases with 
the magnetic lever arm. This is what happens in solutions 
located at the left lower part of Fig. ^ (high lever arms A 
require tiny mass loads k. Fig. ||). 



The other constraint emerges from the requirement 
that jets must become super- Alfvenic. This is expressed 
by ua > 1, which can be written as 

X™(x^-^0>^^^™ (44) 

providing the lower limits in the same figures. These two 
constraints strongly depend on both am and e. 

It turns out that the above two constraints behave dif- 
ferentely with am- As a result, the parameter space shrinks 
considerably (upper and lower limits merging together) for 
decreasing and widens for increasing q,„. We show in 
Fig. |4[ the parameter space for a„i — 0.3; the space is 
considerably reduced and produces no high-w^ jet. Physi- 
cally, it means that if turbulence is not strong enough, the 
magnetic torque will not be able to store enough angular 
momentum in the magnetic field (parameter q) to accel- 
erate matter to the Alfven surface. A much smaller mass 
load would avoid this problem, but this would require an 
overwhelming magnetic pinching and no trans-SM solu- 
tion can be found for cold jets. On the other hand, am 
cannot be much bigger than unity otherwise a too strong 
magnetic compression due to the toroidal magnetic field 
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will occur (Eq. (^2|)). Numerically, we do not found any 
solutions for am > 3. 

No solution with a dominant viscous torque (A <C 1) 
has been found. As showed in Sect. 2.4, super-A jets would 
require magnetic lever arms much bigger than £~^. This 
is impossible here because of the tremendous pinching 
due to the radial component of the magnetic field. Hence, 
both super-SM and super-A constraints merge together 
for decreasing torque ratios A. Thus, cold disc-driven jets 
launched from a large radial extension always carry a large 
fraction of the disc angular momentum (A > 1). Further- 
more, our solutions always produced a large magnetic lever 
arm (A > 3). This allows us to generalize the constraint 
on the maximum ejection index ^ (see Sect. 3.3.2 in F97), 
namely 



Pitch angle 



Opening Angle 



\/l3- 3 A 

4 aTT 



(45) 



which is compatible with the numerical values found. 

Only a rather small range of disc aspect ratios allows 
cold MAES to exist. This can be understood quite easily. 
While a maximum value is imposed by the disc radial 
equilibrium, cold jet production implies a minimum value. 
Indeed, the disc radial equilibrium writes at the midplane 



=S= l-e 



^ + 2(2-/3) 



1/2 



(46) 



where 



= 2gM 



1 + A 



aye{l + A) 



(47) 



is the sonic Mach number. Deviations from Keplerian law 
are mainly due to the radial gradient of plasma pressure 
(of order e^) and the magnetic tension (of order TZm£^)- 
Plasma pressure alone makes it obvious that requiring the 
disc to rotate implies that it cannot be too thick. However, 
the magnetic tension, which is a slowly decreasing function 
of height, provides an even stronger constraint. Indeed, 
as the density falls down, the magnetic support becomes 
more and more important and leads to no rotation (or 
even a negative one) for too thick discs. Numerically, we 
did not found solutions for e > 0.3. 

It is noteworthy that there is no dissipative cold MAES 
for discs too thin (e < 0.05). This is because of the condi- 
tion for launching a cold jet. Indeed, Eq. ( pl| ) implies that 
Xm must be of order s to have A 1. But this cannot 
be verified for very small s because of the trans- Alfvenic 
condition (^J). One can easily see that the thicker disc 
and the more dissipative MAES. 

3.4- Can identical jets be produced from different discs? 

The complexity of Figures I and I raises the question 
whether identical jets could be produced from different 
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Fig. 6. Plots along a magnetic surface anchored at Tq of 
the principal characteristics of two MAES, with the same 
set of jet parameters but different disc parameters (see 
Sect. 3.4). One is a dissipative disc (solid line) where both 
torques are comparable (A = 1.95), whereas the second is 
a weakly dissipative disc (dashed line) with a much higher 
magnetic torque (A = 10.75), hence a larger accretion 
velocity. The pitch angle is defined as arcsm{—B^/ Bp) 
and the opening angle as 9 = a.rccos{Bz/ Bp). In the ideal 
MHD region {z{a) > 0.2ro), the magnetic field configura- 
tion and the velocity field are the same. 



discs. Once in ideal MHD regime, equations (|24| ) to (|28| ) 
and Eq. ^ completely determine the jet structure. The 
leading equation is the Grad-Shafranov equation (p8|), 
which provides the jet transverse equilibrium. In the most 
general case, the solution is elliptic below the Alfven sur- 
face. This means that the sub-Alfvenic solution is com- 
pletely determined once we choose the values of all rele- 
vant quantities at the Alfven surface. In the self-similar 
case, the transfield equation reduces to a second order dif- 
ferential equation of the form 



d'^r Q f dr 
dz"^ T \ ^ dz' 



(48) 



where r{z) describes the shape of magnetic surfaces (Con- 
topoulos & Lovelace 1994 ). Therefore, identical jet struc- 
tures are produced whenever both the leading dimension- 
less parameters {k,X,uja) and conditions at the Alfven 
point are the same. It can be shown (see Appendix B) 
that the locus of the Alfven surface is given by 



To 
ZA 
TA 



cot (jiA — cot ( 



1 - 



U)A 



kA3/2 



cos a A 



(49) 



where za is always found to be of order rA- The magnetic 
configuration (opening and pitch angles) at the Alfven 
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Bo 



Btj, 

Br, 



Ai/2 

K 

UJA 

dr{z) 



dz 



tan Oa{k, \,uja) 



(50) 



(51) 



whereas dynamical quantities verify 



^qTq 

Up 



UJA 



(52) 



wa(i - gA) 



wfiere both the opening angle 9a and the amount qa 
of current available at the Alfven surface {I a /Ism — 
(7A(l+2cr^^) ~ (?a) are only functions of k, A and loa (Ap- 
pendix A). Thus, the three jet parameters impose both the 
location of the Alfven surface and the whole jet behaviour. 

As a consequence, if different sets of disc parameters 
(e, am, 'Pm, Xm) providc the same set of jet parameters 
{k,\,uia), then one gets identical jet configurations from 
different discs. The reason why this occurs lies in the resis- 
tive MHD conditions prevailing inside the disc. Thus, as 
long as there is no theory fixing the parameters for MHD 
turbulence, this situation remains. 

We illustrate this property by showing in Fig. || the 
main characteristics of two jets with roughly (to whithin 
10%) the same jet parameters. The first one (solid line) 
is launched from a dissipative disc with straight mag- 
netic field lines, described by (e, Um, Vm, Xm) — (0.1,1.8,1, 
0.139) and releasing 34 % (A = 1.95) of the total me- 
chanical power. The resulting jets have (k, X,uja) = 
(7 10~^,65,1.56). The second one (dashed hue) is launched 
from a weakly dissipative disc with bent magnetic field 
lines, (e, a^, Vm, Xm) = (0-1, 0.945, 1, 0.366) and releasing 
only 8.5% (A — 10.75) of the mechanical power as radia- 
tion. In this case, the resulting set of jet parameters are 
almost the same, namely {k,\,loa) = (7.8 10"'^, 58, 1.63). 
However, a computation of a "real" MAES requires more 
parameters. Indeed, one must also specify the mass M* of 
the central object, both inner (r^) and external (rg) radii 
where the MAES is established and the accretion rate Mae 
at the external radius. While for most astrophysical ob- 
jects, the central mass and accretion rate can be quite 
severely constrained, the radial extension of a MAES re- 
mains unknown. This is related to the amount of available 
magnetic flux, which must be assumed a priori. Thus, if 
one specifies the same accretion rate, a larger value of A 
(larger magnetic torque) gives rise to a larger accretion ve- 
locity, hence a smaller density at the disc midplane. This 



has several consequences that could be used to discrimi- 
nate between these two solutions: (1) less dense discs can 
be potentially optically thin for the same accretion rate 
(FP95); (2) the ratio of the jet kinetic power to the disc 
luminosity is different; (3) different absolute values of the 
disc magnetic field (/i or Bo) therefore of the jet plasma 
density {pa)- In Fig- 0, we show the cross sections of two 
MAES producing the same jets, namely with same veloc- 
ities, density stratification and magnetic field configura- 
tion. This was accomplished by imposing the same density 
at the disc midplane, which implies a different accretion 
rate. 

To summarize, quasi-identical jet configurations can 
indeed be obtained from different discs. However, physical 
values such as accretion rates (or jet density) and ratio of 
disc luminosity to jet kinetic power, would be different. 
Thus, there is an one to one correspondance between a 
disc and its (cold) jets. 

3.5. Comparison with other steady-state models 

The main assumptions in the work presented here are: (1) 
discs are Keplerian and the launching zone is extended, 
(2) jets are cold. Thus, all models of jets from accretion 
discs obtained with these two basic assumptions should 
verify our analytical constraints. 

Self-similar jets were already obtained by BP82, War- 
die & Konigl ( |1993D and Li ( |1995| ). Basically, our solu- 
tions are compatible with theirs, but we obtained a much 
smaller parameter space. In particular, they were able to 
obtain jets with huge magnetic lever arms A and corre- 
spondingly small mass fluxes k. While BP82 did not treat 



the disc, Wardle & Konigl (1992) replaced mass conserva- 
tion with the prescription pu^ = Constant and Li (1995) 
imposed a static vertical disc structure. A correct treat- 
ment of the disc vertical equilibrium shows that large mag- 
netic lever arms produce an overwhelming magnetic pres- 
sure squeezing the disc and forbidding any steady-state 
solution (see F97 for more details). 

As in BP82, the jets found in this paper always un- 
dergo a recoUimation, until the FM critical point is met. 
This systematic asymptotic behaviour is strongly influ- 
enced by the chosen self-similar geometry (F97). Neverthe- 
less, using exactly the same self-similar ansatz, Contopou- 



los & Lovelace (1994) and more recently Ostriker (1997) 



did found jets with different behaviours. They were able 
to find jets with either an oscillating pattern or reaching 
a cylindrical collimation. This is possible only by vary- 
ing the magnetic flux distribution in the disc, namely for 
/3 > 1. However, disc-driven jets must follow the scaling 
(3 = 3/4 + ^/2 (FP95) where the trans- Alfvenic constraint 
imposes f < 1/2. Thus cold jets with (3 > 1 do not match 
the boundary conditions imposed by Keplerian accretion 
discs. 

To our knowledge, only Li (1996) produced 
magnetically-driven jets from a Keplerian disc with 
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Fig. 7. Poloidal cross-section of two MAES driving jets with the same normahzed jet parameters as in Fig. ^, around 
a one solar mass protostar. Poloidal magnetic field lines are displayed in solid lines, streamlines in dashed lines and 
number density {\ogiQ{nH /cm^^)) is shown in greyscale, from 12 to 7 by 1. In the upper model, the accretion rate, 
chosen to be Mae — W~'^MQ/yr, is provided by both the turbulent viscous torque and the magnetic torque due to 
the jet. As a result, 34% of the accretion power is released as radiation. In the lower model, the magnetic torque due 
to the jet is dominant (only 8.5% of the power is radiated away), hence providing a higher accretion velocity. For the 
same plasma density at the disc midplane as in the upper model, one gets an accretion rate of 3.83 10^"^ AlQ/yr. 



straight magnetic field lines inside. In his model, ambipo- 
lar diffusion is the source of diffusivity in a quasi-neutral 
disc, whereas accretion is solely produced by the magnetic 
torque due to the jet (no turbulent "viscous" torque). 
Thus, almost no power is radiated by the disc and 
Vm << 1- However, his solution displayed fi << 1, 
thereby allowing a high mass loss rate from the disc 
(weak magnetic compression). This is outside our pa- 
rameter range where we found fi to be always of order 
unity. This discrepancy might arise from the different 
set of equations used. In particular, Li ( |1996| ) did not 
use the induction equation for the magnetic field, but 
assumed instead frozen-in ions and got the magnetic field 
behaviour through the prescribed vertical dependency of 
the couphng coefficient. 

4. Summary and concluding remarks 

We investigated the full parameter space of cold, 
non-relativistic, magnetized accretion ejection structures 
(MAES). A generalization of previous works (FP95, F97) 
has been made by taking into account a radial transport 
of angular momentum in the disc. This turbulent trans- 
port is assumed to arise from MHD instabilities triggered 
inside the disc. As a consequence, a cold MAES is char- 
acterized by four local parameters: the disc aspect ratio 
e, the strength of the MHD turbulence am, the magnetic 
Prandtl number Vm and the turbulent dissipation aniso- 
tropy Xm- 

Using self-similar solutions, we were able to obtain 
continuous solutions, from the disc midplane to super- 
Alfvenic jets. The parameter range where cold solutions 



are possible is quite narrow: 10~^ < e < 0.3, 0.3 < am < 
3, < Vm ^ 4 and 10^^ < Xm < 1- Inside this param- 
eter space, both dissipative and weakly-dissipative discs 
are allowed. A dissipative disc releases a significant frac- 
tion of the accretion (mechanical) power as radiation at its 
surfaces. The magnetic torque is comparable to the turbu- 
lent "viscous" torque (A ^ 1). A weakly dissipative disc 
transfers most of the available power to an MHD Poynt- 
ing flux, which is finally converted into jet kinetic power. 
Only a small fraction of the power is radiated away, and 
the magnetic torque is dominant (A ^ 1). 

The importance of the dissipation is controlled by the 
properties of the MHD turbulence. But since a MAES 
must also magnetically drive cold winds, there is a link be- 
tween these properties and jet launching conditions. This 
crucial relation is expressed in Eq. (pT|), namely 

A ^ . (53) 

amav 

Dissipative discs require Xm ~ i-e. anisotropic magnetic 
diffusivities for thin discs. In other words, turbulence must 
lead to a dissipation of the toroidal field stronger than for 
the poloidal field. Dissipative MAES require either am or 
ay bigger than unity. Thus, if one demands those parame- 
ters to be always smaller than unity, only weakly dissipa- 
tive MAES are possible. We believe this is a general result 
for steady-state structures. It arose from analytical consid- 
erations and is thus independent of our self-similar ansatz. 
However, it does depend on two assumptions, namely a 
Keplerian disc and cold ejection. 

Applying these assumptions to an extended launching 
zone, we derived general analytical links between relevant 
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jet parameters {k, A and uja) and disc parameters. The 
most powerful relation, linking accretion to ejection, is 
Eq. @ 



A A3/2 



- Oa) 



1+AA- 1 



sm (pA 



(54) 



which must be bigger than (but of the order of) unity. 
The larger uja, the larger the maximum jet radius. The 
above equation, coming from magnetic flux conservation, 
shows how this jet parameter is intimately related to the 
accretion process {nis = Uo/^Kh). This is quite obvious 
for a dominant magnetic torque, where the toroidal field 
at the disc surface is large. Would the "viscous" torque 
be dominant (A <C 1), the sonic Mach number would be 
very small, TOs ~ a^e. The constraint uja > 1 would then 
require a huge magnetic lever arm from a thin disc. It is 
doubtful such solutions could exist, the disc being prob- 
ably too much squeezed by the magnetic pinching force. 
Another alternative would be a tremendously high viscos- 
ity parameter at, ~ e~^, which is unphysical. 

Two biases were nevertheless introduced in our numer- 
ical solutions. The first one arises from the geometry itself. 
Indeed, self-similar solutions are only valid for MAES set- 
tled on a large radial extension. As a consequence, a cylin- 
drical rather than a spherical geometry is already imposed. 
This fixes the Alfven surface, so that za ~ ta is always 
verified {(j)A between 20° and 50°), hence making grav- 
ity negligible there. As a result, it can be shown that A 
must be larger than (or of the order of) unity. Therefore, 
no cold solution with A <C 1 on a large radial extension 
can be found. However, note that such a small value of A 
would have tremendous implications on the MHD turbu- 
lence parameters Xm, ctm and Uy. 

The second bias was already found by PP92. A class 
of jets described with the same parameters in all magnetic 
surfaces become super-fast magnetosonic (in the conven- 
tional sense, Up > Vfm) and then recoUimate toward the 
axis. Self-similar jets from Keplerian discs fall into this 
category (F97). A realistic modelling of self-coUimated jets 
must therefore be done outside a self-similar framework. 
However, the interrelations between the underlying disc 
and its jets (links between invariants) should be used. 
Ouyed et al. (1997) performed numerical simulations of 
MHD jets from Keplerian accretion discs, using the disc 
as a boundary condition. They obtained solutions that 
recollimate also toward the axis and produce fast MHD 
shocks. Knots are episodically formed at a fixed position 
and, once formed, propagate further down the jet. They 
proposed that this recoUimation and unsteady knot pro- 
duction is realized whenever the toroidal magnetic field is 
large enough, namely when N — B'^/2^oPUp > 1 at the 
disc surface. It is remarkable that our calculations of cold 
MAES satisfy this criterion. 

We showed that there is an univocal link between an 
accretion disc and its jets. But interpreting observations 
is quite tricky, since it depends on assumptions about jet 



emission properties (Cabrit et al. 1999). However, global 
considerations can already provide some constraints on 
theoretical models, like the ratios 2Mj/Ma {— ^Inre/r^) 
and Pjet/Prad (= A). If A is found to be larger than unity, 
then the portion of the disc responsible for mass ejection 
(from ri to r^) should consistently produce a hole in the 
spectral energy distribution. Moreover, since we found an 
ejection efficiency ^ which lies typically around 0.01, mass 
fluxes could provide another measure of the MAES's radial 
extension. Young stellar objects produce jets with typical 
velocities of 300 to 500 km/s. If we assume that these 
jets are produced by a MAES located after the corotation 
radius, too high velocities are obtained. This is because 
the magnetic lever arm is quite large. Thus, models with 
smaller lever arms A should be sought. However, this is 
impossible for cold jets within a self-similar framework. 

Another alternative would be that the emission is dom- 
inated by the external radii of the MAES. In that case, 
only the outer velocities would be inferred from the obser- 
vational data, thus underestimating asymptotic velocities. 
To settle this question, one needs to solve the energy equa- 
tion along the jet and produce synthetic emission maps 
(Garcia et al., in preparation). 

As said previously, we trust our results are quite gen- 
eral for Keplerian accretion discs launching cold jets from 
a large radial extension. The next generalization is to re- 
lax the assumption of cold ejection. Indeed, even if plasma 
pressure may have no dynamical role in jet equilibrium, it 
might produce non-negligible effects on the ejection mech- 
anism itself. In particular, since the vertical equilibrium 
at the disc surface may be deeply modified, different mass 
loads and/or initial jet velocities are expected. The pres- 
ence of hot coronae is quite commonly argued for accretion 
discs, from both theoretical (Galeev et al. 1979| ) and ob- 
servational (e.g. Kwan 1997) grounds. Introducing a hot 
corona is consistent with the presence of a turbulent "vis- 
cous" torque. Indeed, such a torque probably arises from 
3D MHD turbulence in the disc. Numerical simulations 
tend to show the formation of a m agnetized, hot medium 
above the disc (Stone et al. 1996| ). In our model, a small 
fraction of the mechanical power liberated by this torque 
could be converted into coronal heating. This work is cur- 
rently under progress. 

Finally, it has been argued that advection dominated 
accretion flows could also produce self-coUimated jets ( 
see Narayan & Yi 1995, Soria et al. 1997). These are 



quite promising configurations, especially in the context 
of galactic and extragalactic compact objects. In this pa- 
per, we mainly focused our attention on thin (even slim) 
accretion discs. But our framework, with the further de- 
velopment of a hot corona, will allow us to investigate the 
behaviour of flows of that kind. 
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Appendix A: Physical constraints at the Alfven 
point 



If we assume the Alfven surface to be conical, as in radial 
self-similar works or in some numerical simulations of disc- 



In Grad-Shafranov equation (g8|) evaluated at the Alfven 
point (m = 1), the most difhcult term to figure out is the 
density gradient. To solve this problem, we consider the 
function g = 1 — fi/fi*, which can be written 



driven jets (Sakurai 1987, Krasnopolsky et al. 1999), then 
the magnetic flux conservation provides 
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Taking the gradient of this equation and noting that 
Pa/p provides 



where (j)A is the angle between the Alfven surface and the 
vertical axis {za/ta = cot^A in a self-similar solution). 
Injecting this result in Eq. (B.l) gives 
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Thus, the regularity condition at the Alfven point is 
2 
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where we made the assumption that jet invariants are the 
same in every magnetic surface. Inserting the regularity 



condition in the Grad-Shafranov equation (28) provides 
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since J7o — f^* in thin discs. Our solutions always dis- 
played 1 < UJA < 2 and never higher values. The reason 
f or t hat behaviour is hidden in Grad-Shafranov equation 
( A. 5 ), which imposes a highly non- linear link between uja 
and the opening angle 9a- Nevertheless, we observe that 
9a increases with oja (see e.g. Fig. 6 in F97), while the 
above expression of uja shows that it should decrease for 
increasing 9 a- Thus, the jet transverse equilibrium seems 
to impose here some feedback, forbidding high values for 
LOA- As a consequence, Bernoulli equation prohibits small 
values of the magnetic lever arm A. 

Smaller magnetic lever arms could however be ob- 
tained for different geometries. To illustrate this, we ex- 
amine a spherical Alfven surface, in the particular case 
where magnetic field lines are almost straight (Najita & 
Shu 1994). In this case, one has 



This is a general result as long as the jet parameters are 
the same for every magnetic surface. In the case of radial 
self-similarity, it takes the much simpler form 
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where = 2(3 — 3. Thus, the opening angle at the Alfven 
point depends only on gA, k, X and uja- Now, Bernoulli 
equation evaluated at the Alfven surface provides 

^ (A.6) 



A3/2(l + 22/r^)l/2 



2.31 

'- = '--x-:J, 

Note that the last expression is valid for every model of 
disc-driven, cold jets. If the magnetic lever arm is very 
small (A < 3), then the fastness parameter must be large, 
namely loa > 2, even taking into account the gravity term 

(ffi>o). 

Appendix B: Shape of the Alfven surface and uja 



The definition ( p2D of uja can be written, using the defi- 
nition of the mass load k, 



UJA 



VLtf B; 



p,A 



(B.l) 



where i?+ is the poloidal magnetic field at the disc surface. 
Large values of uja with small magnetic lever arms are 
allowed here because of the large dilution of the magnetic 
fiux. 
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